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Abstract 

The GERmanium Detector Array, Gerda, searches for neutrinoless double beta decay 
in 76 Ge using bare high-purity germanium detectors submerged in liquid argon. For 
the calibration of these detectors 7 emitting sources have to be lowered from their 
parking position on top of the cryostat over more than five meters down to the ger- 
manium crystals. With the help of Monte Carlo simulations, the relevant parameters 
of the calibration system were determined. It was found that three 228 Th sources 
with an activity of 20 kBq each at two different vertical positions will be necessary to 
reach sufficient statistics in all detectors in less than four hours of calibration time. 
These sources will contribute to the background of the experiment with a total of 
(1.07 ± 0.04(stat)lo;}g(sys)) x 10" 4 cts/(keV-kg-yr)) when shielded from below with 
6 cm of tantalum in the parking position. 

Keywords: Gerda, neutrinoless double beta decay, calibration, Monte Carlo 
simulation, 228 Th 



1. Introduction 

Experiments with solar, atmospheric and accelerator neutrinos as well as with re- 
actor and accelerator antineutrinos have provided compelling evidence for neutrino 
mixing and for non-zero neutrino masses [see e.g. [1] and references therein]. The ob- 
served phenomena are neutrino oscillations in vacuum and flavour transformation in 
matter. Both phenomena depend on the three neutrino mixing angles, which have been 
experimentally determined [2| [3l H]. Moreover, the former depends, at the first order, 
on the absolute value of the larger squared neutrino mass difference, the latter on the 
smaller one in absolute value and sign [5]. While these quantities are known, the sign 
of the larger mass difference, the absolute value of the masses and the CP violating 
phase (or phases) remain to be measured. 

These experiments are not sensitive to the nature of neutrinos, namely whether neutrino 
and antineutrino are the same (Majorana) or different (Dirac) particles. The observa- 
tion of the neutrinoless double beta decay (Of/3/3) would prove that neutrinos have a 
Majorana component and that consequently lepton number is not conserved. The mea- 
sured half-life will provide information on the Majorana effective electron neutrino mass 
m ee , using the calculated nuclear matrix elements. Several experiments have searched 
for 0u/3/3, but none provided positive uncontroversial evidence [H [H [TUl E] • Cur- 
rent experimental upper limits on m ee are in the range of 0.3 — 1.0 eV. The most 
stringent upper limits on the half- lives have been established by experiments on 76 Ge 
(with enriched HPGe diodes) (Heidelberg-Moscow p2] and IGEX pi]), on 130Te, with 
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natural tellurium oxide bolometers (CUORICINO [H]), on lOOMo with enriched metal 
foils (NEM03 H3HS1) and on 136Xe with a liquid TPC (EXO-200 [H]) and dissolved 
in liquid scintillator (KamLAND-ZEN [18]). 

The GERmanium Detector Array, Gerda, aims to significantly improve these re- 
sults using bare HPGe detectors submerged in liquid argon [19j . To calibrate these 
detectors, radioactive sources are used. This paper presents the Monte Carlo studies 
necessary to determine the optimum type of source, its activity, calibration position and 
background contribution. It is structured as follows: Section [2] gives a brief overview 
of the Gerda experiment, while section [3] specifies the requirements on the calibration 
system. The general setup of the Monte Carlo simulations is described in section [4| 
the results are shown in section [5] In section [6j the background contribution from the 
calibration sources is determined. Section [7] presents the conclusions. 

2. The Gerda Experiment 

The Gerda experiment is searching for the Ou/3/3 decay in 76 Ge which has a half 
life above 10 25 yr. To detect this rare process it is necessary to shield the detectors 
from the radioactivity of the surrounding materials as well as cosmic radiation. For the 
latter the underground location of the Laboratori Nazionali del Gran Sasso (LNGS) of 
INFN was chosen, where 3800 mwe (average for the muon angular distribution and the 
mountain profile) of rock shields the experiment [20]. Furthermore, only materials with 
a very low intrinsic radioactivity were used to build the experiment, and the detectors 
are submerged nakedly into liquid argon (LAr). The argon acts simultaneously as 
cooling liquid and shielding material. If equipped with photomultipliers, it can be used 
as an active veto [21]. A sketch of the cryostat with all relevant dimensions is shown 
in figure [T] The argon cryostat is surrounded by a 10 m diameter water Cerenkov veto, 
which further shields the detectors. 

The experiment is foreseen to proceed in two phases. Phase I is using reprocessed 
HPGe detectors from the Heidelberg-Moscow and IGEX experiments enriched to 86% 
in 76 Ge with a total mass of 17.7 kg together with natural HPGe detectors from the 
Genius Test Facility |22j. The goal is to improve upon current limits on the Qv/3/3 
sensitivity with an exposure of 15 kg- yr and a background of < 10~ 2 cts/(keV-kg-yr)). 
In case no events will be observed above background, a half life limit of Tj/2 > 2.2 x 
10 25 yr can be established, resulting in an upper limit for the effective neutrino mass 
of m ee < 0.23 — 0.39 eV [23J. Phase II will use in addition newly developed broad- 
energy germanium detectors |24j enriched in 76 Ge, and is foreseen to run for a total 
exposure of 100kg-yr. With an aimed background of < 10~ 3 cts/(keV-kg-yr)) a half 
life sensitivity of T]/ 2 > 15 x 10 25 yr and a corresponding effective neutrino mass of 
m ee < 0.09 — 0.15 eV [23] can be reached. 

3. The Calibration System 

For energy as well as pulse shape calibrations of the HPGe detectors radioactive 
sources are used. These sources are staying inside the cryostat to prevent radioactive 
contamination of the liquid argon due to reentry. The parking position of the sources 
during physics runs is on top of the cryostat which is indicated in figure [T] During 
a calibration run the sources are lowered from the top of the cryostat down to the 
detectors. The entire calibration including the movement of the sources should not 
take more than 4 hours to reduce dead time for physical measurements. 

This requires that the activity of the calibration sources should not be too high, such 
as to enable sufficient shielding while in parking position, to ensure a low background 
contribution. On the other hand, the activity has to be high enough to guarantee 
that the total calibration will be performed in a reasonable amount of time. To reach 
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Figure 1: Sketch of the Gerda cryostat. Indicated are the positions of the detector array as well as 
the parking position for the calibration sources. 



this goal Monte Carlo simulations were used to optimize parameters such as activity 
as well as the number of sources and their exact positions during calibration runs. 
Furthermore, the resulting background contribution was determined together with the 
necessary shielding of the source(s). This paper focuses on Phase I of the experiment 
although most of the results can be transferred to Phase II as well. Further details can 
be found in |25|. 



4. Monte Carlo Simulations 

The Monte Carlo simulations were performed in the MaGe |26j framework, a pack- 
age based on Geant4 |27l [28], which simulates the geometry of the entire Gerda 
experiment as well as all relevant physics processes. If not indicated differently, the 
following setup was used: The 8 enriched and 4 non-enriched detectors were placed 
according to the original Phase I configuration shown in figure [2j They are arranged 



in four strings Dl-4 with three detectors per string (see figure 2(a) ). The detectors are 
different in size (see figure 2(b) ) as shown in table [T] they are aligned to the top edge 
of the top detectors. 

Calibration sources at three different positions Sl-3 were simulated with a total of 



10 8 decays; their position relative to the detector array can be seen in figure 2(a) The 
source encapsulation and intended composition were included. If the daughter is not 
stable, the full decay chain was simulated with a pause of 100 microseconds between 
the decay of the different isotopes of the chain to prevent unrealistic pile- up events. 

The resulting energy spectra were folded with the expected energy resolution of the 
detectors which was determined from first test data taken in May 2009 with a 60 Co 
and a 232 Th source. Table [2] shows the energy resolution obtained at different energies 
for one detector. 
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(a) Technical drawing of the lit of the cryo- 
stat showing the position of the four detec- 
tor strings Dl-4 and the three calibration 
sources Sl-3. 



(b) Side view of the detector array. 
Shown in red are the detectors with 
their different heights, shown in grey 
are the three calibration sources with 
their absorbers. 



Figure 2: Configuration of detectors and calibration sources during Gerda Phase I. 



Table 1: Configuration of the detector array including the height h of the different detectors. Dl-4 
refer to four strings with three detectors each. D4.3 is by far the smallest detector and thus is special. 



h [mm] 


Dl 


D2 


D3 


D4 


1 


84 


84 


85 


85 


2 


86 


101 


94 


105 


3 


107 


108 


106 


68 



5. Results 

5. 1 . Type of Source 

Since the Q- value of 76 Ge is at 2039 keV [29], the 5 keV region around this value was 
chosen as the region of interest (ROI). Therefore, several lines in the energy range up 
to 2.5 MeV are necessary for the energy calibration, the requirements are the following: 
There should be at least one line close to the (J-value. A sample of single-site events 
(SSE), which is necessary for pulse shape calibration, can be achieved with a clear 
double escape peak (DEP). This requires a strong full energy peak (FEP) well above 
2 MeV, such that the probability of pair production is sufficiently high, with both 
annihilation photons escaping. A peak to background ratio of at least 2:1 was considered 
sufficient in this case. A full energy peak close to the DEP would be an asset because it 
would allow the comparison of a sample of SSE and a sample of multi-site events (MSE) 
without the influence of the energy dependent energy resolution of the detectors. Since 
the source will stay in the cryostat, a half life of at least several months is required. 

Possible emission of a particles by the source is an important issue because they 
might produce neutrons in (a,n) reactions. The probability of (a,n) reactions depends 
on the energy of the a particles and the threshold energy of the material close to the 
source. Neutrons can contribute to the background in the ROI due to scattering or 
neutron capture, the latter resulting in radioactive isotopes which might emit photons 
or betas with an energy close to the ROI. 

Taking these requirements into account, three possible calibration sources are con- 
sidered: 56 Co, 228 Th and 238 U; table pi summarizes their relevant characteristics. These 



sources were studied in Monte Carlo simulations, positioned as explained in section 5.2 
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Table 2: Energy resolution for different energies taken with a 60 Co and a 228 Th source and one of the 
Gerda Phase I detectors during testing in May 2009. 



E [keV] 


511 
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969 
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1333 


2104 


2615 


a [keV] 
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2.12 


1.78 
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Figure 3: Comparison of different possible calibration sources. Shown is the sum of all entries in all 
detectors normalized to a source activity of 3 x 20 kBq (the motivation for this value is given in sections 
15.21 and 15.31) and the total detector mass. 



Figure [3] shows the sum of the energy spectra of all detectors for the different sources 
normalized to the same activity, time and detector mass. 

As expected, all three sources show several well pronounced lines in the relevant 
energy region. The advantage of 56 Co is that it is not emitting any a particles. Fur- 
thermore, 56 Co has with E 1 = 2.035 MeV a FEP very close but still distinguishable 
from the Q- value. The disadvantages of 56 Co are its short half life of Tim = 77 d and 
no DEP with sufficient statistics for a pulse shape calibration. The next source con- 
sidered is 238 TJ with a half life of T]y 2 = 4 x 10 9 yr and a FEP close to the Q-value 
at E~ = 2.204 MeV. Unfortunately is does no show a DEP with sufficient statistics. 
Moreover it is a possible neutron producer due to spontaneous fission or emitted a 
particles with energies up to E a = 8.8 MeV. The last source considered is 228 Th with a 
half life of T 1( / 2 = 1.9 yr and a peak close to the Q-value at Ely = 2.104 MeV. This is the 
single escape peak coming from the 208 T1 line at E^i = 2.615 MeV. It also results in a 
reasonable DEP at Edep = 1.593 MeV with a peak to background ratio of 2:1. With 
the 212 Bi line at E Bi = 1.621 MeV there is a FEP very close to the DEP which makes 
it ideal for pulse shape calibration. The disadvantage of the source is that it emits a 
particles which can produce neutrons in the surrounding material. The implications 



will be discussed in section 6.2 Nonetheless, it is the best option for the experiment 



and is therefore used for the actual calibration measurements [Tffi [3U] . 
5.2. Positioning 

Scattering processes in the liquid argon lead to a decrease of events in the detectors 
as well as a reduction of the peak to background ratio of the DEP with increasing 
distance between detector and calibration source. Figure [4] shows the total number of 
counts in one detector for different distances to the source. It is evident that the source 
should be as close as possible to the detector. Additionally, the detector strings shield 
each other. Given the string configuration shown in figure [2] it is obvious that more 
than one calibration source is necessary. 

The limiting factor for the positioning of the sources is the lid of the cryostat which 
defines the (horizontal) xy positions of the detector strings and the calibration sources. 
Predefined are the two flanges for the detector strings: One DN250 CF housing three 



detector strings and one DN160 CF housing one detector string (see figure [2(a) ) . The 



smallest usable flange for the calibration sources are DN40 CF giving first limits for 
the smallest possible distance between source and detector. Space is also limited due 
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Table 3: Relevant characteristics for the three possible calibration sources. E max gives the maximum 
energy photons are emitted with. The energy of the line closest to Q = 2.039 MeV with an intensity 
above 1 % is given. It is also indicated if the sources emits any a particles. 





56 Co 


228 Th 


238 tj 




77 d 


1.9 yr 


4 x 10 9 yr 


£ 7 ,max [MeV] 


3.5 


2.6 


2.4 


Ey closest to Q- value [MeV] 


2.035 


2.104 


2.204 


4,max [MeV] 




7.7 


8.8 



l:10' 



20 40 60 80 100 120 140 

d [cm] 



Figure 4: Total number of counts in one detector for different distances to the source. An exponential 
decay was fitted to the data points to determine the mean free path of 2.6 MeV 7's in LAr. The dotted 
line marks the average distance between source and detector as it is realized in Gerda Phase I. 



to the fact that the pipes above the detector flanges have to be accessible to insert the 
detector strings as well as the calibration sources. Furthermore, the lowering system for 
the calibration sources, although reduced in size to an absolute minimum, needs some 
space. Taking all these factors into account, an optimum of three sources positioned 
horizontally as shown in figure [2] is found. The average distance between source and 
detector of 16 cm is marked in figure [4j 

The determination of the best (vertical) z positions for the sources is non-trivial due 
to the different heights of the detectors (see table [I]) , the self-shielding of the detectors 
and an absorber used to shield the sources in their parking position (see section [6]). 
Thus, Monte Carlo simulations were used with the same number of decays for each 
simulation and scanning through the z-space; z = corresponds to the top edge of 
the top detectors. Since every source can be positioned separately they were simulated 
independently. To determine the optimum position for each source, only the relevant 
detector strings were taken into account: For source SI string Dl was considered, for 
S2 the strings D2 and D3 and for source S3 the strings D2 and D4 (for the labels 
see figure 2(a)). The single escape peak (SEP) of 208 T1 was chosen as reference peak 
because with an energy of E = 2.104 MeV it is closest to the Q-value. Figure [5] shows 
the number of counts in the SEP for each detector against the z position. 

Two different cases were considered: 3 different z positions per source, one for each 
detector of a string, and 2 different z positions with the source between two detectors. 
To determine the best positions for both cases, a Gaussian was fitted to the count rates 
in each detector for the different z positions (see figure [5]). The peak positions of the 
Gaussians were used for the 3 position case, the intersection points for the 2 position 
case. Since the sources S2 and S3 have to calibrate two detector strings, the mean 
value between both strings was chosen. The only exception is the bottom detector 
in string D4: Since it is by far the smallest detector its optimum calibration position 
was weighted double when calculating the best position for source S3. The results are 
shown in table HI 

To compare both cases, the total number of counts in the SEP for each detector 



6 



0,1600 



String 1 




"0 50 100 150 200 250 300 350 400 450 

-z [mm] 



0,1200 



(a) String Dl 
String 3 




50 100 150 200 250 300 350 400 450 

-z [mm] 



0,800 
M 

« 700 
o 

" 600 
500 

w 

tj 400 
c 

g 300 
200 
100 

o, 



String 2 


+ Top Det 






+ Middle Det 






* Botto 


n Det 


T 
































A 

































'0 50 100 150 200 250 300 350 400 450 

-z [mm] 



0,1200 



(b) String D2 
String 4 




50 100 150 200 250 300 350 400 450 

-z [mm] 



(c) String D3 



(d) String D4 



Figure 5: Number of counts in the SEP in each detector normalized to a source strength of 20 kBq, 
1800 s calibration time and the detector mass. Gaussian functions wer fitted to determine their peaks 
as positions with maximum count rate (3 calibration positions per source) as well as their intersections 
as optimum positions between two detectors (2 calibration positions). 



Table 4: Optimum vertical position — z in [mm] relative to the top edge of the top detector of the 
calibration sources SI, S2 and S3 (see figure 2(a) I for 2 and 3 vertical calibration positions per source. 

SI S2 S2 SI S2 S3 



50 50 50 



210 215 215 



360 380 365 



3 Positions 



135 130 130 



280 290 320 



2 Positions 



Table 5: Number of counts in SEP with A — 20 kBq per source in [counts/ 1800s] for 3 or 2 calibration 
positions as given in table [4 
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Figure 6: Simulated energy spectrum of the smallest detector in the array after a calibration run with 
3 sources with an activity of 20 kBq each, 2 different z positions and a run time of 30 min per position. 

were calculated using the optimum positions determined above and normalizing to a 
source activity of 20 kBq per source and a calibration time of 1800 s per position. The 
choice of these values will be explained in section |5.3| The goal was to reach at least 
1000 counts in the SEP in each detector. Table [5] shows the results for both cases. The 
different count rates per detector can be explained by their different masses (see figure [5] 
where the normalization was done). Since two calibration positions are sufficient and 
more time efficient, these positions are used for the actual calibration[19| [30] . 

5.3. Activity 

In the next step the minimum source strength necessary to get sufficient statistics 
in all detectors has to be determined. Sufficient is defined in this case as a minimum 
of 1000 counts in the peak as well as a peak to background ratio of 2 : 1. These values 
should be reached using the lowest possible activity within a total calibration run time 
of less than 4 h including the movement of the sources. As a conservative estimate, a 
total moving time of 1 h was assumed for lowering the source 5-6 m from their parking 
position on top of the cryostat down to their two calibration positions and lifting them 
back afterwards. This leaves in total 3h for the calibration itself meaning 1.5 h for 
each position. Since the half life of 228 Th is 1.9 yr and the sources are left inside the 
experimental setup for the full Phase I (1 — 2yr including commissioning) a calibration 
time of 30 min per position at the beginning of Phase I was assumed to ensure sufficient 
statistics also towards the end of this phase. 

To estimate the necessary source activity the results of Monte Carlo simulations 
with each calibration source in its best positions were combined, therefore simulating a 
full calibration run. Again, the SEP was chosen as the reference peak and the number 
of counts as well as the peak to background ratio was determined for each detector. 
Since the third detector in string D4 (D4.3) has by far the smallest count rate, it was 
used reference. 

A simulation of 3 x 10 decays per calibration source and position leads to 950 counts 
in the SEP in detector D4.3 with a peak to background ratio of 2.4:1. Scaling these 
numbers up to 1000 counts was considered as sufficient, leading together with the cali- 
bration time of 1800 s to a minimum activity of 17.6 ±0.6 kBq. Therefore it was decided 
to use sources with an activity of A = 20 kBq for Phase I. Figure [6] shows the energy 
spectrum of the smallest detector corresponding to this configuration. 

6. Background Contribution 

Since a low background environment is crucial for the success of the experiment, 
an estimate of the background contribution in the ROI from the calibration system 
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in the parking position is essential. Several isotopes in the 228 Th decay chain decay 
via a decay; hence it was considered both, the emitted 7's and the neutrons produced 
via (a,n) reactions in the surrounding material. The background contribution from 
neutrons during physics runs as well as calibration runs will be determined. 

6.1. Gamma Background 

According to original planing, there are a minimum of 3.4 m of liquid argon between 
the sources and the top detectors; the rest is gas and will not be considered here. In 
case of the 7 background, two methods were used and compared. The first one is a 
combination of an analytical estimate with Monte Carlo simulations (MCS), the sec- 
ond one is a pure MCS with the sources in their parking position. The latter is very 
CPU consuming. Therefore the first method was used for faster results to be able to 
determine the necessary shielding. The full MCS were used to verify the results. 

In the semi-analytical method, MCS are used for the last 1 m of LAr. The flux at 
this position is estimated using linear attenuation <p = <J)q x e~ d l l with <^>o being the 
initial flux, d the thickness of the absorbing material, I = ^- the mean free path, [i the 
mass attenuation coefficient and p = 1.394 g/cm 3 the density of liquid argon. The 2.6 
MeV line from 208 T1 is the only line above the Q-value and therefore the only possi- 
ble source for background in the ROI. In the following, only photons with this energy 
will be considered. The NIST database provides experimentally determined values for 
the mass attenuation coefficient for 2.044 MeV and 3.0 MeV photons; the value for 
the 2.6 MeV photons is interpolated according to their fit function. A mass attenua- 
tion coefficient of /i = 0.0359 cm 2 /g was found |31j . resulting in a mean free path of 
^nist = 20.0 ± 0.4 cm. 

To confirm this result within the Gerda geometry, MCS were used: A beam of 
10 s 7's with E = 2.6 MeV was directed to one detector with varying distances between 
them. Fitting an exponential decay to the total counts in the detector over distance 
lead to a mean free path of /mcs = 20.65 ± 0.05 cm. The corresponding plot can be 
found in figure |4j With this method, the mean free path could be determined with 
much higher precision than from the NIST database. Furthermore, the used geometry 
reflects the actual situation in the cryostat and it is more conservative. Since it is also 
more conservative, the mean free path determined by MCS will be used in the following. 

For the initial flux, the number of decays from 3 calibration sources with A = 20 kBq 
in one year was used and then taken into account that due to the branching ratio just 
36% of these decays end with an emission of a 2.6 MeV 7. Since the source radiates 
isotropically but the detector array with a radius of r^et = 15 cm covers only a small 
area of this sphere, this flux reduces by a factor of ^ S phere/^det = 2055, resulting in 
(j)Q = 3.3 x 10 8 7/yr. After 2.4 m of LAr, the flux reduces to <j> = 3.0 x 10 3 7/yr. 

To obtain a conservative estimate of the background in the ROI, a photon beam 



directed to the center of one of the detector strings was simulated. Figure 2(a) shows 
the xy position of the detector strings and the calibration sources. SI pointed to 
Dl, S2 to D2 and S3 to D4. In total, 4.8 x 10 9 7's with an energy of 2.6 MeV 
were simulated. These simulations result in a total of 8.4 x 10 4 events in the ROI. 
Rescaling to the initial flux and using a total mass of 17.7 kg, a background index of 
(2.9±0.1(stat) ±0.3(sys)) x 10" 4 cts/(keV-kg-yr) was obtained. Such a level is tolerable 
for Phase I with a background goal of 10~ 2 cts/(keV-kg-yr). However, since other parts 
of the experiment, especially the cryostat, will contribute to the background as well, 
further shielding is preferred and will be discussed later. 

In the full MCS three calibration sources were placed 3.4 m above the detector array. 
To acquire more statistics in a shorter time, more detector strings were included in the 
simulations leading to a total of 48 detectors with a total mass of 76.8 kg. Again, only 
2.6 MeV 7's are simulated since they are the only possible source for background in the 
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Figure 7: Energy deposition of 2.6 MeV 7's started 3.4 m above the detector array. A total of 1.05 x 10 12 
7's were simulated in a 2n hemisphere which corresponds to the number of 228 Th decays in about 3 years. 



Table 6: Screening results for tantalum 



Isotope 


238tj 


232 Th 60 Co 40 K 137 Cs 182 Ta 


Activity [mBq/kg] 


< 11 


<9 < 1.9 <33 < 2.5 52 ±5 



ROI. Photons emitted in the top 2tt hemisphere will loose too much energy on their 
way to the detectors if they reach it at all and can be ignored. Therefore, just photons 
emitted in the lower 2tt hemisphere were simulated. The simulation of 1.05 x 10 12 7's 
lead to a total 19 events in a 400 keV ROI from 1839-2239 keV. This large region was 
necessary due to the low statistics. The energy spectrum in the full energy range as 
well as the ROI is shown in figure [7J This corresponds to a background contribution 
of (2.0±0.6(stat) ±0.2(sys)) x 10~ 4 cts/(keV-kg-yr). As expected, the semi-analytical 
approach shows the more conservative limit but both values agree within errors. 

To further shield the sources, the possible shielding material needs to be as radio- 
pure as possible with very good absorption properties and has to be machinable. The 
material fulfilling these requirements best is tantalum with /i(2.6MeV) = 0.04cm 2 /g, 
p = 16.7 g/cm 3 and a natural radioactivity of about 50 mBq/kg from 182 Ta with a half 
life of 114 d. The detailed screening results measured with the Gator screening facil- 
ity [32] can be found in table [6j 

A reduction of the 7 background to B~ f < 10 -5 cts/(keV-kg-yr)) was considered 
necessary. Using again linear attenuation to determine the needed minimum thickness 
of the shielding material, d > 5 cm was found. It was therefore decided to shield the 
calibration sources with 6 cm of tantalum to be conservative, leading to a background 
contribution due to the 7 radiation of the calibration sources of 

5 7 = (4.8 ± O.l(stat) ± 0.3(sys)) x 10~ 6 counts/ (keV-kg-yr) (1) 

6.2. Neutron background 

Neutrons from (a,n) reactions in the calibration source can contribute to the back- 
ground in two different ways: By (n,n'7) inelastic scattering of fast neutrons in the 
surrounding materials, especially in argon and due to neutron capture resulting in ra- 
dioactive isotopes emitting photons with an energy close to the ROI. Elastic scattering 
is in this case not important since the maximum energy which can be transferred is far 
below the Q- value of 2039 keV. For the scattering, the contribution from the parking 
position is most important. For the case of neutron capture, the calibration runs are 
also important to consider, since isotopes close or in the detectors might be activated. 
Therefore Monte Carlo simulations for both cases were performed. 

The neutron spectrum was determined using SOURCES4A |33j with the following 
assumptions: A Th02 solution is coated on a gold foil placed in a stainless steel cap- 
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Figure 8: Predicted neutron spectrum from (a,n) reactions of a 228 Th source embedded in gold. 



sule [34J. Due to the high (a,n) threshold of gold, the only reaction partners for the 
a's from 228 Th is oxygen which was taken into account in natural isotopic abundance 
in the calculations. The resulting neutron spectrum is shown in figure [8j 

As a first step, the background contribution from elastic scattering of neutrons 
in the detectors as well as (n,n')7's during the calibration run was inspected. As 
figure [9] shows, this background has no significant influence on the calibration spectrum. 
Therefore it will be ignored in the following. 

For the estimation of the background due to neutron scattering, Monte Carlo sim- 
ulations were started with the sources in their parking position and the neutron spec- 
trum shown in figure [8j The simulation of 10 9 neutrons resulted in 230 counts in 
the energy region between [1950,2060] keV. As a realistic estimate, a neutron flux of 
A n = 10- 3 n/(s-kBq) [M] was assumed. Using these numbers for three calibration 
sources with an activity of A = 20 kBq each the following background contribution was 
found: 

^scatter = ^ 2 ± 0-4 ( stat )) x 10 -5 ± 2 0%(sys) cts/(keV-kg-yr) (2) 

which is sufficiently low for Gerda Phase I. 

In the next step the isotopes produced due to neutron capture reactions in all 
parts of the experiments were analyzed. Figure [TO] shows the type and amount of 



produced isotopes in one year during calibration runs (figure 10(a) ) as well as from the 



parking position (figure 10(b)), assuming one calibration run of lh per week. Most of 
these isotopes do not emit photons nor /3's with energies above 2 MeV and are thus 
no potential background candidates. Produced isotopes emitting photons or /3's above 
2 MeV are: 19 0, 28 Al, 41 Ar, 42 K, 49 Ca, 55 Cr, 56 Mn, 66 Cu and 77 Ge. 

Concerning the photon emitters, those isotope, whose intensity is below 1 %, were ig- 
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(a) Energy spectrum during a cal- (b) Background contribution (c) Background contribution 

ibration run: The 7 contribution from the sources in parking po- from the sources in parking 

(black) is shown together with the sition due to neutron scattering position due to 77 Ge produced 

background from elastric scatter- events. via neutron capture, 
ing of neutrons in the detectors as 
well as (n,n')7's (red). 

Figure 9: Neutron background during calibration (a) and from parking position (b) and (c). 
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(a) Activated isotopes during calibration. 



(b) Activated isotopes from source in parking posi- 
tion. 



Figure 10: Activated isotopes normalized to amount produced per year. 



nored. This reduces the possible background candidates to 9 Ca, 56 Mn and 77 Ge. 56 Mn 
is produced in several volumes, the highest production rates are the rock of the labo- 
ratory (6.2 x 10 4 isotopes/yr) and the wall of the cryostat (2 x 10 4 isotopes/yr), assuming 
again a neutron flux of A n = 10 -3 n/(s-kBq) and a total activity of 
3 x 20 kBq. Due to the continuous neutron flux from the sources in parking position and 
the short half-life of 56 Mn of 2.6 h, the production and decay rate are in equilibrium, 
leading to an activity in the rock of 2.8 x 10 -4 Bq and in the wall of the cryostat of 
8.9 x 1CT 5 Bq. A total of 1.2 x 10 9 decays were simulated in the cryostat wall. No events 
were found above 1.5 MeV leading to an upper limit on the background contribution of 
< 1.5 x 10~ 9 cts/(keV-kg-yr), which is negligible. Since the water tank further shields 
the detectors from the radiation of the rock, its contribution can be ignored. The same 
argument holds for 49 Ca produced also in the rock of the laboratory with a very low 
production rate of 26 isotopes/yr. 

Due to the short stopping range of /3's, only isotopes produced inside the cryostat 
were considered. This reduces the possible background candidates to 41 Ar and 77 Ge. 

Ar, with a half life of 109m and a Q- value for the /3-decay of 2.5 MeV, is produced 
in the liquid argon with a production rate of 2.6 x 10 6 isotopes/yr, assuming again 
a neutron flux of A n = 10~ 3 n/(s-kBq) and a total activity of 3 x 20 kBq. A total 
of 2.4 x 10 9 decays were simulated, isotropically distributed in the cryo liquid. To 
achieve reasonable statistics, the 200 keV interval around the Q- value was chosen as 
ROI. A total of 48 counts were found, leading to a negligible background contribution 
of (7.3 ± 1.1) x 10- 8 cts/(keV-kg-yr). 

Therefore, only 77 Ge, produced in the detectors themselves, is potentially danger- 
ous. After the neutron capture, the 77 Ge is highly excited and de-excites via a 7 cascade 
into either the ground state or the metastable state 77m Ge. The ground state j3 decays 
with T1/2 = 11.3 h, with endpoint energies up to 2.5 MeV and emits several gammas 
with energies up to 2.3 MeV. 77m Ge decays with T1/2 = 53 s with a 19% chance into 
the ground state emitting a 160 keV 7 or with a 81 % chance via (3 decay with endpoint 
energies up to 2.9 MeV, emitting several gammas with energies up to 1.7 MeV. The 
total production rate of both, 77 Ge and 77m Ge, is An Ge = 8.6 isotopes/yr, assumed 
to be equivalent to the decay rate. Since the corresponding cross sections vary signifi- 
cantly EH! EZ] , the range of the background contribution was estimated using two 
simulations: One assuming that the total amount of the produced 77 Ge will decay from 
the ground state and the other one from the metastable state. In both cases a total 
of 2 x 10 7 decays were simulated, resulting in -B gr0 und = 2.9 x 10~ 5 cts/(keV-kg-yr)and 
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-Bmeta = 4.0 x 10 5 cts/(keV4cg-yr). As a conservative limit, £? gr ound wm be used: 

B?7 Gc = (4.0 ± O.l(stat) l?i(sys) ) x 10~ 5 counts/ (keV-kg-y) (3) 

Combining the different background contributions due to neutrons, 

B n = (1.0 ± O.l(stat) ^(sys)) x 10~ 4 counts/ (keV-kg-y) (4) 

was found, which is well below the Gerda Phase I goal but relatively close to the Phase 
II goal. Thus, further reduction is preferred for this later stage of the experiment. 

7. Conclusion 

The Monte Carlo studies for the calibration system of the Gerda experiment 
showed that three 228 Th sources with an activity of 20 kBq each are necessary to cal- 
ibrate the detector array in Phase I of the experiment. To prevent scattering of the 
7's in the liquid argon which is used to cool and shield the detectors, the calibration 
sources were placed in the horizontal plane as close as possible to the detectors. In 
the vertical direction two positions between the detector layers are necessary to reach 
sufficient statistics in each detector within a calibration time of 30 min. 

During a physics run the calibration sources are parked on top of the cryostat 
and the radiation of the sources might contribute to the background in the region 
of interest. Both gamma as well as neutron radiation as result of (a,n) reactions 
were considered and a total background contribution of (1.07 ± 0.04(stat)^Q'Jg(sys)) x 
10~ 4 cts/(keV-kg-yr) was found. This is well below the background goal for Phase I of 
10~ 2 cts/(keV-kg-yr) but might become more relevant for Phase II with a background 
goal of 10 -3 cts/(keV-kg-yr). The highest contribution is due to (a,n) neutrons, which 
make 93% of the total expected background. A shielding of each source with 6 cm of 
tantalum was necessary to reach such a low background contribution. 
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